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Vortex Tiling in a Spin-2 Spinor Bose-Einstein Condensate 
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We investigate core structures of fractional vortices in a spin-2 spinor Bose-Einstein condensate 
and find that commensurability between discrete symmetries of the order parameters inside and 
outside of a vortex core brings about spontaneous breaking of the rotational symmetry of the core 
and the associated twisting excitations. Our predictions can be tested experimentally by using a 
rotating spin-2 87 Rb Bose-Einstein condensate and enable an experimental realization of the cyclic 
order parameter, regardless of whether the ground state is cyclic or otherwise. 
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Quantized vortices are topological defects of the su- 
pcrfluid order parameter. For single component Bose- 
Einstein condensates (BECs), the quantized vortices 
emerge as phase singularities, and the order parameter 
vanishes at the core center. However, the situation dras- 
tically changes for spinor BECs, where the vortex core 
can be filled with a state having different symmetry. For 
example, a spin-1 polar BEC supports a half-quantized 
vortex with a ferromagnetic core and a spin-1 ferro- 
magnetic BEC can hold a spin vortex with a polar core. 
The vortex core structure for the spin-1 BEC is, however, 
uniquely determined since there are only two possible 
ground-state phases in zero magnetic field. The prob- 
lem becomes much more involved in higher-spin BECs 
because of the greater possibilities of the ground-state 
phases. Since the order parameter of the spinor BEC is 
generally anisotropic, the following questions arise which 
we refer to as "vortex tiling problems": (1) which state 
fills the vortex core?; (2) which orientation the order pa- 
rameter of the vortex core takes with respect to the sur- 
rounding BEC?; (3) how can the vortex core be embed- 
ded in the BEC if they have different discrete symme- 
tries? We will show that the commensurability in spin 
space gives rise to breaking of the rotational symmetry 
of the vortex in real space, thus allowing us to probe 
the symmetry of the order parameter by observing the 
deformed shape of the vortex core. 

In this Letter, we address the tiling problems for frac- 
tional vortices in a spin-2 BEC. The ground state of the 
spin-2 BEC is shown to have four different phases under 
zero magnetic field @, H, 0, d, [3| ; ferromagnetic, uniax- 
ial nematic, biaxial nematic, and cyclic, where the last 
three phases have discrete symmetries. The biaxial ne- 
matic phase has the non-Abelian fourth dihedral (-D4) 
symmetry jf|, and the cyclic phase has the non-Abelian 
tetrahedral symmetry [7[. We consider here two highly 
nontrivial cases; a vortex in the biaxial nematic phase 
with the cyclic core, and a vortex in the cyclic phase 
with the biaxial nematic core. The solutions of these 
tiling problems are shown in Fig. [T] Figures Q] (a) and 
(b) show the cloverleaf and triad profiles of the biax- 
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FIG. 1: (Color) (a) and (b) Plots of |$(0,0)| 2 for the bi- 
axial nematic phase and the cyclic phase $c, where 
$(0, (f>) = J2m=-2 Yim{6, 4>)ip m - Colors represent the phase of 
$(0,0). (c) and (d) Distributions of $(0,0) for the 1/2-1/4 
vortex with the cyclic core and the 1/3-1/3 vortex with the 
biaxial nematic core. Here 1/2 and 1/4 refer to the ir gauge 
transformation and the ir/2 spin rotation around the vortex, 
respectively, (e) and (f) Profiles of | ^4_3o | 2 and IA20I 2 for the 
1/2-1/4 and 1/3-1/3 vortices, where the size of each panel is 
16 x 16 in units of the healing length £ = h~/\/2Mc2n to t- 



ial nematic and cyclic order parameters in spin space. 
The direction of the triad (cloverleaf) inside the vortex 
must be determined so as to be smoothly connected with 
the cloverleaf (triad) of the BECs outside of the vortex; 
Fig. Q] (c) (Fig. Q] (d)) shows the tiling pattern of the 
order parameters for these vortices. Each lobe of the 
triad (cloverleaf) in the core is connected with that of 
the cloverleaf (triad) outside the core. The smoothly 
connected directions are indicated by arrows with col- 
ors of the corresponding lobes. As a consequence, the 
rotational symmetry of the vortex core is spontaneously 
broken and the core is deformed to become triangular 
(elliptic) for the vortex in the biaxial nematic (cyclic) 
phase with the cyclic (biaxial nematic) core, as shown in 
Fig. ffl(e) (Fig. ffl(f)). Here the three-fold and two-fold 
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symmetries of the cores reflect those of the cyclic and bi- 
axial nematic phases, respectively. Details of Fig. Q] are 
explained below. 

Another important finding of this Letter is that the 
cyclic state can be realized as the vortex core state in 
an F = 2 87 Rb BEC @, 0, H E|- Although the 
ground state of this BEC has yet to be determined be- 
tween the cyclic and biaxial nematic phases due to the 
complications arising from the quadratic Zeeman effect 
and hyperfine-spin-exchanging relaxations, we have con- 
firmed that the vortex with the cyclic core can be gener- 
ated by rotating the biaxial nematic initial configuration 
of the F = 2 87 Rb BEC, regardless of which phase the 
true ground state is. The core deformation predicted in 
this Letter can also be used as a hallmark of the cyclic 
state. 

We consider a BEC of spin-2 atoms with mass M whose 
mean- field energy functional is given by 2, 3] 
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where ip m is the order parameter of the BEC in a 
magnetic sublevel m — 0, ±1,±2 at position r, and 

n tot = X^m=— 2 IVVrel ) F — ^2im,m' = ~2^mfm,m''4 ! m', 

and A 2 o = (2ip2ip-2 — 2i/>iV'-i + V'oVVS are the total 
number density, the spin density, and the spin-singlet 
pair amplitude, respectively, with f m ^ m > being a vec- 
tor of spin-2 matrices. The ground state of Hamilto- 
nian ([T]) is (i) ferromagnetic (F 7^ and A20 = 0) 
for c\ < and C2 > 20ci, (ii) uniaxial nematic or bi- 
axial nematic (F = and A20 ^ 0) for C2 < and 
c 2 < 20ci, and (iii) cyclic [F = A 20 = 0) for ci > 
and C2 > 0. A representative order parameter for each 
phase is given by (i) $ F = ^/n to t(l,0, 0, 0, 0) T , (ii) *u = 
7^(0, 0, 1, 0, 0) T or B = y/fh^(l, 0, 0, 0, l) T /\/2, and 
(iii) * c = \A 7 ^t(0,\/2,0,0,l) T /v / 3. In the mean-field 
theory, \1/b and Wtj are degenerate under zero magnetic 
field. It is predicted that zero point fluctuations lift this 
degeneracy [6]. The quadratic Zeeman effect also lifts 
this degeneracy, stabilizing * B for the F = 2 87 Rb BEC 
at a magnetic field above the order of mG, which is con- 
sistent with experiments d, 0, [l^, 11 1 . We therefore do 
not consider the uniaxial nematic phase here. 

The order parameter can be expressed in terms of 
the spherical harmonics: <!>(#, 4>) = ^ = _ 2 ^ m (^)^"" 
where 9 and <f> are the polar and azimuthal angles in 
spin space, respectively. Figures Q] (a) and (b) show 
${0, <f>) for the biaxial nematic phase and the cyclic 
phase V^c- The profiles of <&(#, <fi) for the biaxial ne- 
matic and cyclic phases are cloverleaf and triad, reflecting 
the Di and tetrahedral symmetries, respectively. Be- 
sides the trivial 2ir gauge transformation and the 27T 
rotation in spin space, $(0,0) remains invariant under 



some discrete spin-gauge transformations. In the biax- 
ial nematic phase, for example, the combined operations 
of the 7r gauge transformation and the it/2 spin rota- 
tion around the z axis do not change <!>(#, </>), namely 
e iT e -tf~T 7 2 \]/ B = v]/ B . This discrete symmetry admits a 
fractional vortex around which the global gauge changes 
from to 7r and the cloverleaf rotates by tt/2. We refer 
to such a vortex as the 1/2-1/4 vortex, where 1/2 and 
1/4 refer to the tt gauge transformation and the tt/2 spin 
rotation around the vortex, respectively. Similarly, we 
can define the 1/3-1/3 vortex in the cyclic phase [1,0]. 
The order parameters around straight 1/2-1/4 and 1/3- 
1 /3 vortices on the z axis can be expressed in cylindrical 
coordinates (r, ip, z) as 



*i/ 2 -i/4(r -» 00) = eWe-"'*'**! 



*l/3-l/3(r ~> OO) 



= ^(l,0,0,0,e^) T , 

e i<p/3 e -if z <p/3-q, c 



^(0,V2,0,0, e ^f, 



(2) 
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where the vortex is placed at r = 0. Among 
the vortices with nonzero mass circulation k = 
h/M§dl £ m Im[V4(VVm)]/n t ot, the 1/2-1/4 and 1/3- 
1/3 vortices are the most stable ones in the biaxial ne- 
matic and cyclic phases respectively Q- 
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FIG. 2: Schematic diagram of core structures for 1/2-1/4 
(left) and 1/3-1/3 (right) vortices on the ci - 62. plane. F: 
ferromagnetic core, B: biaxial nematic core, and C: cyclic core. 
The range of (ci, c 2 ) for a 87 Rb BEC is indicated [ll[. 

Let us now examine the core structure. The core struc- 
ture is determined by the competition among four terms 
in Eq. (TTJ and strongly depends on the interaction pa- 
rameters co,i,2- We here consider the case for the 1/3- 
1/3 vortex in the cyclic phase: (i) When ci^ en, the 
third and fourth terms in Eq. (fT]) play a dominant role 
in determining the core structure. Because finite F or 
A20 drastically increases the energy, n tot vanishes at the 
core, (ii) When *C Co and c\ < C2, on the contrary, 
ntot is constant with the fixed total number of particles 
N = J dx n to t because of the inequality / dx n% ot ^ N 2 . 
A nonzero F raises energy less than A20 because of 
Ci < C2, and the core is ferromagnetic, (iii) When 
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ci.2 Co and ci > C2, the core is biaxial nematic for non- 
vanishing A 2 o with constant n to t- (iv) When a.2 ~ Co, 
the last three parts in Eq. |T]) compete. To minimize 
the kinetic energy, components with no phase singularity 
must be constant. Actually, m= 2, 1,0,-1 m Eq. ^ are 
constant inside and outside the core and the order param- 
eter becomes v / 2n tot /3(0, 1,0,0,0) at the ferromagnetic 
core. Similarly, we can show that the core of the 1/2-1/4 
vortex is cyclic for | cx,a | "C Co and c\ > 0. Figure [2] shows 
a schematic diagram of the core structure in the ci - 62 
plane, where c~i = Ci/cq and 62 = c^jc^. 

We next discuss the orientation of the order parameter 
of the 1/2-1/4 vortex with the cyclic core that has to be 
embedded in the biaxial nematic order parameter. The 
energetically most favorable configuration of the order 
parameter at the core is found to be 

*i/a-i/ 4 (r = 0) = e 2 */ 3 e-^/ 3 e-^c 

= Vntot/3(l, 0, 0, V2e iC , 0) T , (4) 

where £ is an arbitrary real number. Other possible 
candidates such as e^l^P cos ~ 1{1/ ^ ) ^i/ 2 -i/i{r = 
0) = (i/2, 0, l/\/2, 0, i/2) are energetically unfavor- 
able, because it would require a drastic change from 
#1/2-1/40" -> 00) to *i /a _i/4(r = 0). In Fig. [T] (c), 
we show $(0,0) for the 1/2-1/4 vortex with the cyclic 
core. As shown in the figure, Eq. (J4|) fixes the order 
parameter of the core state in a manner such that the 
(0, 0, Indirections in Figs. [1] (a) and (b) coincide with 
each other. As £ changes, the order parameter shown 
in Fig. [1] (b) rotates around the z axis accompanied by 
a gauge transformation. Since the m = — 1 component 
vanishes at r — > 00, £ does not contribute to the vortex 
energy. The remaining SO(2) symmetry of the core state 
corresponding to £ is spontaneously broken, leading to a 
localized Goldstone mode in the vortex core. For the 
1/3-1/3 vortex with the biaxial nematic core, we obtain 
the entire $(0,0) distribution as shown in Fig. Q] (d). 
In this case, the (1, 1, 0)-direction in Fig. Q] (a) and the 
(0, 0, Indirection in Fig. Q] (b) coincide with each other. 
We now consider the ansatz for the 1/2-1/4 vortex 



the transverse magnetization is given by 

■(ae^, 5 ,/3e^, 7 ,/e^) T , 
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(6) 



2 and 
where = 

1+6^2^-2) is 



resulting in the two-fold anisotropy, where a, /3, 7 = 0, 
/ = 1, and g = y2 at r ~ > 00, and a = /3 = / = 
0, g = VW2, and 7 = y/ij2e K at r = 0. We have 
numerically minimized the mean-field energy ([1]) with the 
ansatz $1/2-1/4 and $1/3-1/3 with £ = 0. Figures Q] (e) 
and (f) show the obtained profiles of |^3o| 2 and 1^20 
for $1/2-1/4 and $1/3-1/3, respectively, 

-9(V> 2 ^-2 + V'-i^s) + + &M>- 

the spin-singlet trio amplitude which is a measure of the 
cyclic order Q. We can clearly see the triangular and 
elliptic deformations of the cores. 

In Fig. Q] (c), the triad of the cyclic order parameter 
inside the core is embedded with the cloverleaf of the 
biaxial nematic order parameter, and there are three di- 
rections along which one of the three lobes of the triad is 
smoothly connected with one of two lobes of the clover- 
leaves surrounding the core. Along these directions, the 
triad is well overlapped with cloverleaves, which corre- 
spond to the minimum of |_F| 2 and the vertices of the 
triangle profile of | 1 2 - Similarly, the biaxial nematic 
core of the 1/3-1/3 vortex is deformed to be elliptic as 
shown in Fig. [T] (d). Thus, the triangular and elliptic 
shapes of the vortex core are the manifestations of the 
discrete three-fold and two-fold symmetries of the triad 
and the cloverleaf. Since the order parameter of the fer- 
romagnetic phase does not have any discrete symmetry, 
vortices with the ferromagnetic core are isotropic. 
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where a,/3, 7 = and / = g = 1 at r — > 00, and 
a = (3 = / = 0, and g = ^/2~Jl and 7 = 2/%/3e< at r = 0. 
Phase singularities for m = 1,0 components are intro- 
duced to minimize the nonzero transverse magnetization 
The obtained transverse magnetization becomes: 
F 2 +F 2 = n tot |2( 5 a*+ 7 /*)+V6(/37*+«/3*e 3 ^)| 2 /6, and 
the three-fold anisotropy arises. With this anisotropy, 
the Gross-Pitaevskii equation [1] derived from Eq. 
also becomes anisotropic, leading to the anisotropic den- 
sity distribution of each magnetic sublevel |0 m | 2 . Simi- 
larly, the ansatz for the 1/3-1/3 vortex which minimizes 



FIG. 3: (Color) (a) |y4 30 | 2 and (b) |0 O | 2 of the 1/2-1/4 vortex 

and (d) 1 02 1 2 of 



with the coupling constant c% , and (c) \A2o 
the 1/3-1/3 vortex with cf\ (e) \A S 

(ii) 



with cf>. (f) \A 



2 • v c ; 1^-30 
2 of the 1/3-1/3 vortex with dtp 



2 of the 1/2-1/4 vortex 
Here a^o 



is the harmonic oscillator length. \A3o\ and other values are 
normalized by n t ot(r = 0) 3 / 2 and ntatir = 0), respectively. 

We show that the vortices shown in Fig. [T]can be gen- 
erated in the F = 2 87 Rb BEC under an external rota- 
tion. For the F — 2 87 Rb BEC, ci ~ (8.80 ±0.53) x 10~ 3 
and c~2 — (— 4.71±5.15) x 10~ 3 were experimentally ob- 
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tained [xi| . from which we can anticipate the 1/2-1/4 
vortex with the cyclic core or the 1/3-1/3 vortex with 
the biaxial nematic core in this system (see Fig. 
They have nonzero mass circulations k = h/2M and 
h/3M, respectively, and can be experimentally gener- 
ated by an external rotation. In the presence of a trap- 
ping potential V = M[uj 2 (x 2 + y 2 ) + lj 2 z 2 ]/2 and an 
external rotation iMl z (xd y — yd x ) along the z axis, we 
find the stationary state of the order parameter for the 
F = 2 87 Rb spinor BEC by minimizing the total en- 
ergy. With experimental values of the radial frequency 
uj r and the two-dimensional density [l(| , we perform two- 
dimensional simulation under Q z = 0.15w r . We use Co i 2 
measured in [llj ]. For C2, we adopt the following two 
cases; c 2 (i) = -4.71 x 1CT 3 and c 2 (ii) = 0.44 x 1CT 3 
where the ground state is expected to be biaxial nematic 
and cyclic, respectively. The stationary state can be ob- 
tained by calculating the imaginary-time development of 
the Gross-Pitaevskii equation derived from Eq. ([1]) 
We choose two initial states; stationary <3>b and \&c with- 
out rotation for the 1/2-1/4 vortex and the 1/3-1/3 vor- 
tex, respectively. The results are shown in Figs. [3J (a) 
- (d). Both 1/2-1/4 and 1/3-1/3 vortices are stabilized 
under the rotation, and the spontaneous breaking of ro- 
tational symmetry also occurs reflecting the symmetry 
of the core state, which is shown in Fig. [3] (a) for |^3o| 2 
in the 1/2-1/4 vortex and Fig. [3J (c) for \A 20 \ 2 in the 
1/3—1/3 vortex. This symmetry breaking brings about 
an anisotropic density distribution, in particular, for the 
component described by a and /3 in Eqs. © and (J6|), and 
our results can be experimentally tested by the Stern- 
Gerlach experiment. In the present simulation, \ipo\ m 
the 1/2-1/4 vortex and \ip2\ 2 in the 1/3-1/3 vortex are 
highly anisotropic as shown in Figs. [3] (b) and (d). 

We have numerically confirmed that the vortex tiling 
phenomenon occurs independently of the sign of C2. Even 
when we start from the "false" ground state like \I/b with 
or with c| , those states involve the "metastable" 
1/2-1/4 or 1/3-1/3 vortex, because the values of c\,i are 
so small that the spin relaxation dynamics from the false 
to true ground state is much slower than that of the vor- 
tex nucleation. Figures [3] (e) and (f) show the 1/2-1/4 
vortex in 'J'b with c 2 n ' and the 1/3-1/3 vortex in \fc 
with c 2 , and those figures are quite similar to Figs. [3] 
(a) and (c), respectively. These results imply that the 
nucleated vortex depends only on the initial state regard- 
less of the sign of c 2 , and both cases of the 1/2-1/4 and 
1 /3— 1 /3 vortices can be experimentally realized using the 
F = 2 87 Rb spinor BEC. In particular, the cyclic order 
parameter can be realized as the core of the 1/2-1/4 vor- 
tex surrounded by the biaxial nematic order parameter. 
Since both 1 /2-1 /4 and 1/3-1/3 vortices are non-Abelian 
[H> 0] > the above method can be used to produce the first 
experimental realization of non-Abelian vortices. 



The vortex tiling leads to a new type of excitation of 
vortices, that is a "twisting mode". Due to the Gold- 
stone mode with respect to ( in Eq. ([4]), the triangular 
and elliptic vortex cores can rotate without energy cost. 
When we consider the 3-dimensional vortex line, there 
should be the characteristic wave excitation with a lin- 
ear dispersion which corresponds to the twisting of the 
triangle or elliptic shape along the vortex line and can 
be defined as "twisting wave" . The detailed study of this 
mode will be reported elsewhere. 

In conclusion, we have studied fractional vortices in a 
spin-2 spinor BEC and shown that the problem of vor- 
tex tiling in the order parameter space arises. A 1/2-1/4 
vortex (1/3-1/3 vortex) with a cyclic core (biaxial ne- 
matic core) in the biaxial nematic (cyclic) phase requires 
the vortex tiling between the discrete -D4 symmetry of 
the biaxial nematic state and tetrahedral symmetry of 
the cyclic state. The rotational symmetry of the cores is 
spontaneously broken and their shapes become triangular 
for the 1 /2-1 /4 vortex and elliptic for the 1 /3-1 /3 vortex, 
which directly relate to the three-fold and two-fold rota- 
tional symmetries of the triad and cloverleaf structures. 
The vortex tiling can be experimentally implemented by 
using a rotating F — 2 87 Rb BEC. 
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